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^T ■ The conductance resulting from resonant tunneling through a droplet of 

N ~ 30 electrons is used to measure its chemical potential /i/v- Abrupt shifts 
of ujv occur at sharply defined values of the magnetic field, at which the 

state of the droplet changes. These are used to study part of the phase- 

O ■ 

diagram of the droplet in strong magnetic fields; we find evidence for a new 



cd 



S^ . phase in the spin polarized regime. We make a detailed comparison between 

theory and experiment: Hartree-Fock provides a quantitative description of 

the measurements when both spin-split states of the lowest orbital Landau 

level are occupied and a qualitative one in the spin polarized regime. 
PACS numbers: 73.20.Dx, 73.20.Mf 
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Recent experiments |l],|2| have demonstrated the possibility of measuring the chemical 
potential /xjy of a droplet of N electrons confined by an external potential, an artificial atom. 
Abrupt shifts of {In occur at values of the magnetic field B at which the ground state (GS) of 
the droplet changes. These results have stimulated calculations of the B-N phase diagram, 
in which each phase is designated by the quantum numbers of the GS; the changes in nn{B) 
happen at the phase boundaries. Because exact numerical calculations are possible only 
for N < 6 IHH], approximate methods [g|-0] have been used for larger iV to account for 
electron-electron interactions. The strong magnetic field regime is the appealing place to 
test these approximations, because the most intriguing aspects of the phase-diagram occur 
at these fields. In particular, MacDonald et al. || , and Chamon et al. have independently 
predicted the existence of new phases of a spin polarized droplet in a parabolic potential. 
These phases are especially interesting because any transition in the spin polarized regime 
is the consequence of many-body phenomena that cannot be explained by a single-electron 
picture. 

In this Letter, we present detailed measurements of a portion of the phase diagram in 
strong magnetic field. We propose a new systematic approach for comparing the experi- 
mental results with the models. Using this approach, we find that Hartree-Fock (HF) [|5],[7J 
provides a quantitative description when both spin states of the lowest orbital Landau level 
(LL) are occupied, whereas a semi-classical model (SC) || does not, indicating that ex- 
change plays an important role. We also find evidence for a new phase in the spin polarized 
regime, which is described qualitatively by HF. 

The device that we study is of the type described by Meirav et al. ||. It consists of 
a two-dimensional electron gas (2DEG) in an inverted GaAs/ALGai_ x As heterostructure 
with electrostatic gates above and below it. The bottom gate is a highly conducting sub- 
strate of n + doped GaAs. A positive bias, V g , applied to the bottom gate varies the density 
of the 2DEG. On the top surface of undoped GaAs, two metallic (TiAu) gates are litho- 
graphically patterned with a double constriction. Applying a negative bias to these top 
gates depletes the 2DEG lOOnm underneath them, confining the electrons to an island be- 



tween the constrictions. Current flows through the resulting electron droplet via the tunnel 
barriers caused by the constrictions. The top gate geometry of the device under investi- 
gation has been examined with an Atomic Force Microscope. The constrictions are poorly 
defined, but we estimate that the region between them is roughly 500 x 500nm 2 . According 
to the simulation of the device by Kumar et al. [[J, the external confinement potential of 
the droplet is approximately parabolic. Although all results presented here are for this one 
structure, we have observed similar features in samples of different geometries. 

The negative bias on the top gate is maintained constant during the experiment and 
the bottom gate voltage is varied in a narrow range near V g = 160 ± lmV, for which 
the electron density of the 2DEG regions outside the constriction is almost constant at 
1.3±0.01xl0 11 cm -2 . The conductance G of the device as a function of V g at B = 0T is shown 
in the lower inset of Fig.l. It consists of quasi-periodic sharp peaks {AV g = 1.2mV), crudely 
described by the coulomb blockade mechanism [|T0| | . In this model, when the bottom gate 



voltage is set between peaks, transport is suppressed by the charging energy U ~ 0.66meV 
necessary to add an electron to the droplet. Each period thus corresponds to the addition 
of one electron to the droplet. At resonance, the electrochemical potential of the droplet, 
fiN — eaV g , is aligned with the Fermi energy of the leads and current flows; current requires 
a fluctuation of the charge on the droplet. Thus, the value of V g at which the peak occurs 
provides a measure of fijy. At T= 0, fiN = E N — E N _i, where E N is the energy of the 
TV-electron GS. 

We begin by considering the effect of magnetic field on a single conductance peak. That 
is, we measure ^n{B) at constant N . The value of the gate voltage at which the iVth 
conductance peak occurs is plotted as a function of magnetic field between 1 and 5T in 
Fig.l. McEuen et al. recognized that the change in behavior near 1.6T results from the 
depopulation of all but the lowest orbital LL ||. The step-like behavior of the peak position 
above 1.6T can be thought of as resulting from the transfer of electrons between the two spin- 
split states of the lowest orbital LL [J5J. Each step corresponds to a change in the quantum 
numbers of the GS, for example, the total spin of the droplet. The number of steps above 



1.6T is proportional to N, the number of electrons in the droplet, but the proportionality 
constant depends on the shape of the charge distribution. 

A novel way to characterize the data in Fig.l is to examine the separation in B of the 
upward steps. The peak conductance as a function of B has a sharp minimum at each of 
these steps [fIJ]. The nth minimum precisely determines B n , the field for the nth step. (The 
B n are indicated by arrows in Fig.l.) We plot in Fig.2a the quantity (B n — i? n _i) _1 as 
a function of B n . Because each step corresponds to the flip of a spin, one may think of 
(B n — Bn-i)^ 1 as being roughly proportional to the spin susceptibility. A fit to the form 
y {B) = y [(B - B')/B'Y gives B' = 1.7±0.02T and e = -0.41 ±0.06 for our data; the solid 
curve in Fig. 2a shows the fit. The same functional form also fits the experimental data for 
two other devices with different geometries (500 x 700nm 2 and 450 x 900nm 2 ) and larger 
numbers of steps (~ 25 and 35); we find e = —0.37 ±0.1 for all three devices. 

Plotted in Fig. 2c is the result obtained when (B n — i?„_i) _1 is determined using /!jv of 
the SC model of McEuen et al. [[J. For a parabolic potential with cylindrical symmetry, 
V(r) = m*c<j 2 r 2 /2, the SC spatial density of electrons p(r) is approximately that of classical 



electrons in zero magnetic field p(0Wl — (r/R) 2 , except near r corresponding to integer 
filling factors v = 2ir£ 2 p(r) where the electrons form an incompressible liquid, (m* is the 
effective mass of electrons in GaAs, uoq is the oscillator frequency, R is the radius of the 
droplet, and £ is the magnetic length.) In this picture, the change in behavior at B = 1.6T 
corresponds to a filling factor v = 2 at r = 0. At fixed B, p{r) is uniquely determined by 
N and Tiojq. We adjust iV so that the calculated uat(-B) has the same number of steps as 
observed experimentally. With N fixed, we adjust Tiuj® to match the value of B' = 1.7T, at 
which the transfer of electrons between the two spin states of the lowest LL begins in our 
experiment. Using this procedure, we find N = 42 and Jiujq = 1.8meV. As seen in Fig. 2c the 
SC model predicts values of (B n — Bn^i)^ 1 which are roughly the same size as the measured 
ones. However, it does not predict the upward curvature of (B n — 5 n _i) _1 near B' . In the 
SC calculation, the last spin flip occurs at Bj = 3.23T= 1.95'. In fact, the SC value of the 
ratio Bi/B' is almost independent of iV and Tiujq. Although we find a step at 2B' in our 



measurement, there is also an additional step at B c = 3.75T= 2.2B' (filled circle in Fig.l) 
not predicted by the SC model. For all devices studied we find a step at 2B', marking the 
complete depopulation of the higher energy spin state, and a step at larger field, in the spin 
polarized regime JTT]. For one device, we have explored the phase diagram beyond 2.7 B' 



and have found evidence for other steps JT2 

The step at B c behaves in a way that is very different from those between B' and 2B' '. 
By examining successive peaks in G vs. V g (lower inset of Fig.l), i.e. probing the droplet at 
successive N, we find that each step in /i^ (Fig.l) shifts to higher B when another electron 
is added to the droplet. We have averaged the shift over four consecutive conductance peaks 
and have plotted its inverse [B n (N) — B n (N — 1)] _1 in Fig.3a for each of the steps in Fig.l. 
This quantity measures the slope of the phase boundary (dN/dB). It is clear from Fig. 3a, 
that the step at B c has a larger slope than those at lower B. 

The temperature dependence of the step at B c is also peculiar. Figs. 4a and b show that 
the features between 1.7 and 3.4T disappear by 500mK as T is increased. This behavior is 
now well understood JT|. In clear contrast, the height of the step at B c = 3.75T does not 
change with temperature up to 800mK, our measurement limit. 

To compare our results with a more sophisticated theory, we have performed a HF 
calculation of /j,n(B), choosing the states of the symmetric gauge as the complete basis 
set, with the Hilbert space truncated to the two spin states of the lowest LL. Because of 
exchange, the HF p(r) is more compact than the SC p(r) with larger incompressible regions, 
smaller compressible regions and a more rapid decrease with r near the edge of the droplet 
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In Fig. 2b, we plot (B n — 5 n _i) _1 as a function B n extracted from pn{B) for the HF 
model with N = 27 electrons and hu>o = 2.1meV. As with the SC model, these parameters 
are chosen to match the number of steps in the peak position and the experimental value of 
B', respectively. Like the value of Tiojq used for the SC model (1.8meV), the hu>o that fits the 
HF calculation is in agreement with the value 2 ± lmeV calculated by Kumar et al. from 
the sample geometry ||. The HF value of iV = 27, however, is different from that (N = 42) 



which fits the SC model. This discrepancy is a result of the difference in shape of the charge 
distribution. In the HF model, the number of steps in ix^{B) is equal to N/2 because the 
two spin states are equally occupied at B', and half the electrons flip their spin as the field 
is increased between B' and B\. 

It is obvious from Fig. 2b that the HF calculation is in excellent quantitative agreement 
with the experiment (Fig. 2a). This is particularly impressive since there are no other fitting 
parameters once N and huo are fixed. In particular, the HF model predicts correctly the 
apparent divergence of (B n — 5 n _i) _1 near B', in clear contrast with the SC model. A fit 
to the HF results with y(B) gives e = —0.43 ± 0.03, which is the same as the experimental 
value within the errors. The apparent divergence of (B n — Bn^i)^ 1 in Fig. 2a suggests that 
because of exchange the two spin states of the lowest LL are equally occupied at B' in our 
droplet. This is consistent with another experimental observation: a new step in /j,^(B) is 
added between B' and 2B' for every two electrons added, implying that the two spin states 
of the lowest LL of our droplet are equally populated with increasing N. 

The HF calculation predicts that the last spin flip occurs at Bj = 3.15T= 1.855'. Like 
the SC model, the HF ratio Bj/B' is nearly independent of iV or hu . Thus, above Bj the 
droplet is spin polarized. MacDonald et al. 0, and Chamon et al. [|7| showed that there 
exists a region in the B-N phase diagram (sketched in the upper inset of Fig.l) in which, 
for N < N c ~ 100, the GS of the spin polarized droplet is the maximum density droplet 
(MDD). In the MDD state, all the single-particle eigenstates of angular momentum index 
m = 0, 1, ..N — 1 are occupied, leading to an approximately constant p(r) in the droplet. 
The MDD is of course the GS of non-interacting electrons at high B, but surprisingly it is 
also the GS in a region of B-N even in the presence of repulsive interactions |||7]]. 

With increasing magnetic field the radius of the MDD decreases, the electrons get closer 
together, and the interaction energy eventually favors a larger area droplet. HF J/J predicts 
that, at B c , the edge undergoes a reconstruction and electrons form an annulus at a distance 
~ 2£ away from the central droplet, causing an abrupt upward shift of (x^ at B c of roughly 
the same height as the step at Bj J7||. In the HF calculation, B c /Bi is almost independent 



of huo, but it decreases with increasing iV ||[| for N < N c . 

The excellent quantitative agreement between HF and the experiment for B < 2B', 
strongly suggests that the MDD is formed in our experiment above 2B'. The HF calculation 
predicts that the transition to the reconstructed droplet occurs at B c = 4.21T for our droplet, 
a value larger than the one observed experimentally. In this regard, it is important to bear 
in mind that although the HF energy of the MDD is exact because the MDD is an exact 
eigenstate of the many-body Hamiltonian ||, the HF energy of the reconstructed droplet 
is only variational. Therefore, the calculated value of B c is an upper bound on the true 
transition field. Indeed, an exact calculation for small N H|J shows that the HF model 
overestimates B c . 

Turning to the slopes of the phase boundaries, one sees in Fig. 3 that [B n (N) — B n (N — 
1)] _1 from HF (~ 2.2 x 10 _3 G _1 ) agrees fairly well with experiment (3 ± 1 x 10 _3 G _1 ) 
between B' and 2B' . However, at B c , the HF value 3.2 x 10 _3 G _1 , is smaller than the 
experimental value, 8 ± 1.5 x 10 _3 G _1 . The quantities [B n (N) — B n (N — 1)] _1 at Bj and B c 
are the slopes of the phase boundaries in the B-N phase diagram between which the MDD 
is the GS. The fact that [B n (N) — B n (N — 1)] _1 is larger at B c than at Bj suggests that 
the MDD does not exist above some N c ||[7| . The experimental observation of both a larger 
value of [B n (N) — B n (N — 1)] _1 at B c and a smaller value of B c than the ones predicted by 
HF suggests that N c is smaller than predicted by HF. 

We have extended the HF calculation to obtain excited states and thus study the tem- 
perature dependence of /i^. We find that the HF excitation spectrum (proportional to the 
height of the /itv steps) has an energy scale 4 times larger than the experimental one over 
the entire magnetic field range. Nonetheless, the steps in the region between B' and Bj are 
predicted to wash out more rapidly with increasing T than the one at B c in agreement with 
observation (Fig. 4). 

Finally, we note that HF also describes the B dependence of the conductance peak height 
12]] . HF predicts the experimentally observed ||T2[ decrease in peak height just below B c 



followed by an increase for B above B c . The increase above B c is ascribed to the reduced 



separation between the edge of the droplet and the leads when the annulus is formed. 

The failure of HF to predict the size of the magnetic field window in which the MDD is 
the GS (Fig. 2b) and the dependence of B c on N (Fig.3b) may indicate that correlations are 
playing an important role in this transition. The downward step at about 3.5T (Fig.l) is 
also reminiscent of features predicted to result from correlations || . 

In conclusion, we have made a detailed study of the conductance peak positions in strong 
magnetic fields. We have focused on that part of the phase diagram in which only the lowest 
orbital LL with its two spin-split states are occupied. By looking at the increase in magnetic 
field required to flip each successive spin, we are able to make a quantitative comparison 
between experiment and theory. We find that HF is in excellent quantitative agreement 
with experiment at low field. However, when the droplet is spin polarized a new transition 
occurs which is only qualitatively described by HF. 
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FIGURES 
FIG. 1. Upper Inset: B-N phase diagram of the droplet. The boundaries corresponding to a 

change of the total spin of the droplet in the 2 > v > 1 regime are omitted; the MDD domain 

of stability is limited on one side by B[(N), the boundary of the spin polarized phase, and on 

the other side by B C (N), where there is a reconstruction of the charge density. Above N c the 

MDD phase is terminated. Lower Inset: Conductance through the island as a function of the 

bottom gate voltage at B = OT. Main: Position of the Nth conductance peak as a function of B 

at T = lOOmK. We have used a constant factor a = 0.55 to convert the bottom gate voltage scale 

to energy [1,5]. The arrows indicate the minima of the conductance peak height. B n is the field 

for the nth minimum above 1.6T, with n = {1, ..., 14}. 

FIG. 2. (a) (B n — i? n _i) -1 vs. the B n obtained from Fig.l. The error bars represent the spread 
of the data when the analysis is repeated for other conductance peaks on the same device, (b) and 
(c) results obtained with the HF (N = 27 and huj = 2.1meV) and SC (N = 42 and hu> = 1.8meV) 
calculations. The solid line is a fit with y(B) = y [(B — B')/B'] e where B' = 1.7 ± 0.02T and 
e = —0.41 ± 0.06 for the experiment and e = —0.43 ± 0.03 for the HF. Bj indicates the field onset 
of the spin polarized regime in both models. The solid circle indicates B c . The dashed line in (c) 
is the constant interaction model [5] (the scale is expanded by a factor 10). 

FIG. 3. (a) [B n (N) — B n (N — 1)] _1 is the slope of the phase boundary; it is measured by looking 
at the shift in 5-field of the same B n between two adjacent conductance peaks (or equivalently 
when an electron is added to the droplet). The plotted value is the average over 4 consecutive peaks 
and the error bars are the standard deviations (b) HF and SC values of [B n (N) — B n (N — 1)] _1 
measured from the simulated hn(B). 

FIG. 4. Magnetic field dependence of the peak position at 100 and 500mK. (a) and (b) show 
the behavior below 2B', while (c) and (d) show it above 2B' . In each figure, the peak position is 
offset for clarity. 
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